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We report on a search for gravitational waves from coalescing compact binary systems in the 
Milky Way and the Magellanic Clouds. The analysis uses data taken by two of the three LIGO 
interferometers during the first LIGO science run and illustrates a method of setting upper limits 
on inspiral event rates using interferometer data. The analysis pipeline is described with particular 
attention to data selection and coincidence between the two interferometers. We establish an obser- 
vational upper limit of 7?. < 1.7 x 10^ per year per Milky Way Equivalent Galaxy (MWEG), with 
90% confidence, on the coalescence rate of binary systems in which each component has a mass in 
the range 1-3 Mq. 

PACS numbers: 95.85.Sz, 04.80.Nn, 07.05.Kf, 97.80.-d 
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I. INTRODUCTION 

The Laser Interferometer Gravitational- Wave Obser- 
vatory (LIGO) is an ambitious US initiative to detect 
gravitational waves from astrophysical sources such as 
coalescing neutron stars and black holes, spinning neu- 
tron stars, and supernovac. The LIGO detectors arc laser 
interferometers with light propagating between large sus- 
pended mirrors in two perpendicular arms. They mea- 
sure the strain (differential fractional change in arm 
lengths) produced by gravitational waves from astrophys- 
ical sources by monitoring the relative optical phase be- 
tween light paths in each arm Q . LIGO comprises three 
detectors housed at two geographically distinct locations: 
in Hanford, WA, there are two interferometers, one with 
arms 4 km long (which is referred to as HI in this arti- 
cle) and one with arms 2 km long (112); in Livingston, LA 
there is one interferometer with arms 4 km long (LI). The 
LIGO interferometers 0,0] form part of a worldwide net- 
work of gravitational-wave detectors which includes the 
British-German GEO 600 detector the French-Italian 
VIRGO detector 0, the Japanese TAMA300 detector 
and five resonant-bar detectors . 

Among the most likely sources of gravitational waves 
accessible to earth-based detectors are binary systems 
containing neutron stars and/or black holes When 
they reach design sensitivity, the initial interferometers 
in LIGO should be sensitive to gravitational waves gener- 
ated during the last several minutes prior to coalescence. 
Current wisdom suggests that binary neutron star co- 
alescences could provide up to ~ 1/4 events per year 
detectable by the initial LIGO interferometers at design 
sensitivity 0| . Binary black hole coalescences could 
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provide up to ~ 2 events per year The rates, however, 
are uncertain and may be significantly lower. 

Previous published searches for gravitational waves 
from compact binaries used data from the LIGO 40m 
prototype ^ and early data from the TAMA300 detec- 
tor 0. The 40m data was taken in 1994 over a week- long 
run which yielded 25 hours of data and resulted in an up- 
per limit rate of 0.5 events per hour in the Galaxy. The 
instrument was sensitive to sources up to 25 kpc away 
with signal-to-noisc ratio equal to 10. The TAMA300 
data was taken in 1999 over three nights which yielded 
6 hours of data and resulted in an upper limit of 0.59 
events per hour for events producing a signal-to-noise ra- 
tio larger than 7.2, corresponding to sources up to 6 kpc 
away. Searches for generic gravitational- wave bursts have 
also been performed using data from multiple detectors 
which operated simultaneously. Over 100 hours of data 
froni prototype interferometers at Glasgow and Garch- 
ing |l2l |. and four years of data from the International 
Gravitational Event Collaboration (IGEC) of resonant- 
bar detectors resulted in event rates consistent with the 
background of the instrumental noise 0, . 

This article reports on the first search for gravitational 
waves from binary neutron star inspiral using LIGO data. 
The first scientific data run, called SI, lasted 17 days in 
2002 and involved all three LIGO detectors. The detec- 
tors were sensitive to binary inspiral events to maximum 
distances (at signal-to-noise 8 in a single detector) be- 
tween 30 and 180 kpc, depending on the instrument, al- 
lowing the most sensitive search yet. (The TAMA300 
collaboration is currently analyzing ~ 1000 hours of 
data which will provide a comparable upper limit.) The 
GEO 600 detector collected data in coincidence with 
LIGO during the entire SI run and achieved an excellent 
duty cycle of 98%. At the time of SI, GEO 600 was still 
being commissioned and was operated without signal re- 
cycling - an essential part of its final optical design. It 
was therefore operating at a sensitivity significantly lower 
than that of the LIGO detectors and its own target sensi- 
tivity. Hence GEO 600 was not included in this analysis. 
The upper limit reported here, TZ < 1.7 x 10^ per year 
per Milky Way Equivalent Galaxy (MWEG), is the best 
direct observational limit on binary neutron-star coales- 
cence to date. This rate is far from expected astrophysi- 
cal rates, but demonstrates the progress of instrumental 
commissioning and success of the data analysis effort. 

Many of the analysis techniques presented here will 
be used in future searches for gravitational waves. For 
instance, we expect to use these methods while analyzing 
data taken during the second LIGO science run between 
February and April 2003 when the detectors had roughly 
ten times better amplitude sensitivity than in SI. 

The outline of the paper is as follows. Section |H] con- 
tains a description of the instruments, performance, sen- 
sitivity and duty cycle during SI. Section lTTll describes in 
detail the target population of binary neutron-star sys- 
tems and the gravitational waves they generate. The 
matched filtering technique used to search for these sig- 
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nals in the data is reviewed in Sec. IIVI Filter outputs 
above a certain signal-to-noise ratio threshold constitute 
triggers which are cataloged for further analysis, provided 
they satisfy a test to determine the consistency of 
the data with the expected waveform. Section |3 de- 
scribes data quality cuts and instrumental vetoes which 
are applied to eliminate triggers from times when the rel- 
evant interferometer was not operating properly. Surviv- 
ing triggers are passed through an analysis pipeline which 
generates a list of event candidates from a combination 
of multi- and single-interferometer data, as detailed in 
Sec. IVII To avoid statistical bias, the veto conditions 
and pipeline parameters were tuned using a playground 
data set which was representative of, but separate from, 
the main data set. An upper limit on the rate of binary 
neutron star coalescences is calculated in Sec. IVIII and 
systematic errors are considered in Sec. IVIIll Section Hxl 
summarizes the results and discusses the prospects for 
future data runs. 



II. THE LIGO DETECTORS 

The LIGO interferometer design is a variant of a 
Michelson interferometer, with a laser light source and 
a beam splitter which directs the light along two perpen- 
dicular arms. Mirrors at the ends of the arms reflect the 
light beams back to the beam splitter, where they re- 
combine and interfere according to their relative optical 
phase; this interference provides a sensitive measure of 
the length difference between the two arms. To augment 
the basic Michelson design, partially transmitting input 
mirrors are placed near the beam splitter to form a long 
Fabry-Perot cavity in each arm with a finesse of ~ 220. 
An additional partially transmitting mirror is placed in 
the path of the input laser beam to form a composite 
power-recycling cavity, which increases the amount of 
light circulating in the interferometer. A more detailed 
description of the LIGO optical configuration and other 
instrumentation may be found in Ref. |14| . 

The light source for each interferometer is a medium 
power Nd:YAG laser, operating at a wavelength of 1.06 
^m . Before the light is directed into the interferome- 
ter, its frequency, amplitude and direction are stabilized 
using a combination of active and passive stabilization 
techniques. 

To isolate the mirrors and other elements from ground 
and acoustic vibrations, the detectors employ active and 
passive seismic isolation systems |l6l IT^ , from which the 
mirrors are suspended as pendulums. These form a cou- 
pled oscillator system with high isolation for frequencies 
above 40 Hz. The mirrors, major optical components, 
vibration isolation systems, and main optical paths are 
all enclosed in a high vacuum system. 

Various feedback control systems are used to keep the 
multiple optical cavities tightly on resonance and well 
aligned [T^l. The strain signal s{t) = [Lx{t) — Ly(t)]/L is 
derived from the error signal of the feedback loop used to 



control the differential motion of the interferometer arms. 
To calibrate the error signal, the effect of the feedback 
loop gain is measured and divided out, and the response 
R{f) to a differential arm strain is measured and factored 
in. The absolute scale of the response is established us- 
ing the laser wavelength by measuring the mirror drive 
signal required to move through a given fraction of a 
fringe. The response varied over the course of the SI run 
due to drifts in the alignment of the optical elements; it 
was tracked by injecting fixed-amplitude sinusoidal sig- 
nals (calibration lines) into the differential arm control 
loop, and monitoring the amplitudes of these signals at 
the measurement (error) point po| . 

The interferometer noise is characterized by the one- 
sided power spectral density Sn{f) of the signal s{t). The 
sources of noise that are expected to limit the eventual 
sensitivity of the LIGO detectors are shot noise (deter- 
mined by circulating light power, dominant at high fre- 
quencies), thermal noise (determined by energy dissipa- 
tion mechanisms in the mirrors and suspensions, dom- 
inant at intermediate frequencies), and seismic noise 
(dominant at low frequencies). Figure ^ shows the ex- 
pected noise due to these effects (at LIGO's design tar- 
get), expressed as RMS strain noise, along with typical 
spectra achieved by the LIGO interferometers during the 
SI run. (Typical GEO 600 noise during SI is also shown 
for comparison.) The differences among the three LIGO 
spectra reflect differences in the operating parameters 
and hardware implementations of the three instruments 
which are in various stages of reaching the final design 
configuration. For example, all interferometers operated 
during SI at a substantially lower effective laser power 
level than the eventual level of 6 W at the interferometer 
input. Thus the shot-noise region of the spectrum, above 
200 Hz, is much higher than the design goal. In addition, 
the SI configuration only had a partial implementation 
of the laser frequency and amplitude stabilization sys- 
tems, and a partial implementation of alignment control 
systems for the mirrors and the beam splitters. Despite 
these shortcomings, the detectors were sensitive to bi- 
nary neutron star coalescences within the Galaxy and 
the Magellanic Clouds as illustrated in Fig. |21 

The 17-day run yielded 363 hours of data when at least 
one interferometer was in stable operation. The three 
interferometers were simultaneously in stable operation 
for 96 hours. For the analysis presented in this article, 
we chose to use data only from the two 4 km detectors, 
LI and HI. While H2 was nearly as sensitive as HI, 
its noise exhibited a greater degree of non-stationarity, 
leading to a rate of spurious triggers which would have 
compromised the sensitivity of the search. LI and HI 
were simultaneously operational for 116 hours during the 
SI run, providing data for the first combined analysis 
of interferometric detectors sensitive to inspiral events 
throughout the Galaxy. In addition, they were separately 
operational for 54 and 119 hours, respectively. 
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FIG. 1: Typical sensitivities of tlie LIGO a nd GEO 600 interferometers during tfie SI data run, sfiown as equivalent RMS 
strain amplitude spectral density /!.rms(/) ~ \J fSn{f), where Sn{f) is the one-sided noise power spectral density. Typical noise 
spectra for the two 4km interferometers, LI and HI, used in our analysis are shown in the left panel; the smooth solid curve 
indicates the target sensitivity of the LIGO 4 km interferometer design. Spectra for the 2km interferometer H2 and GEO 600 
are shown in the right panel; the smooth solid and dashed curves indicates the target sensitivities of the LIGO 2 km and 
GEO 600 interferometer designs. The thick lines with arrowheads show the characteristic strains, hchai{f) ~ fh{f), expected 
from binary neutron star systems (optimally located and oriented with respect to the detector) during the last few minutes 
before coalescence. These characteristic strains are approximately equal to the amplitude of a gravitational wave signal at a 
given frequency times the square-root of the number of cycles produced in a logarithmic band about the given frequency. The 
ratio of /ichar to /irms in the sensitive band of the instrument provides an estimate of the signal-to-noise ratio that could be 
achieved in detecting such a signal using matched filtering. When the LIGO instruments are operating at the target sensitivity, 
inspirals of double neutron stars (2 x 1.4Mq) are expected to be detectable within an equivalent volume « (47r/3) x (21 Mpc)^. 



III. TARGET POPULATION AND 
WAVEFORMS 

Radio observations of pulsars confirm the existence of 
binary neutron star systems in the Galaxy |2ll . . Gen- 
eral relativity predicts the decay of a binary orbit due to 
the emission of gravitational radiation. The decay rate 
inferred from observations of PSR1913-I-16 agrees with 
the prediction within 0.3% [Hilllil. The orbital de- 
cay is easily modeled for compact binary systems con- 
taining neutron stars or stellar mass black holes. The 
binary orbit is expected to evolve through the LIGO fre- 
quency band by the emission of gravitational waves alone, 
making it possible to accurately compute the evolution 
without reference to complicated micro-physics. 

When a compact binary system first forms, the or- 
bit may be widely separated and highly eccentric. (See 
Rcf. Q for a discussion and plots of birth separations 
and eccentricities). Gravitational radiation, emitted pre- 
dominantly at twice the orbital frequency of the binary 
system, causes the orbit to shrink and circularize (much 
faster than it shrinks p^) so that the binary components 
eventually spiral together along a sequence of nearly cir- 
cular orbits with decreasing period. For binary neutron 
stars or stellar-mass black holes, the gravitational ra- 
diation eventually enters the frequency band of earth- 
based gravitational-wave detectors. At this point, the or- 



bit decays rapidly and the gravitational waveform chirps 
upward in frequency and amplitude, sweeping through 
LIGO's sensitive band. During SI, the LIGO interfer- 
ometers were sensitive to gravitational-wave frequencies 
above about 100 Hz; an inspiral signal from two IAMq 
objects would traverse the sensitive band in 2 seconds. At 
design sensitivity, the sensitive band will stretch down to 
~ 40 Hz and the signals will spend about 30 seconds in 
the sensitive band. 

For low-mass binary systems, the waveforms are well 
approximated by a post-Newtonian expansion [2^l2Ml29l | 
in the LIGO frequency band. Due to the uneven con- 
vergence of this expansion and a still indeterminate co- 
efficient at higher order, we used second-order post- 
Newtonian waveforms |2a | in this analysis. The wave- 
forms are parameterized by the masses of the two com- 
panions / = (mi, TO2), the inclination of the orbit relative 
to the plane of the sky, and the starting orbital phase. 
Other orbital parameters such as eccentricity and spin 
are not expected to be significant for binary neutron star 
coalescence |^ (s^, so we do not consider them in 
this analysis. The strain produced in the instrument is 
written as 

h{t) = ^-^^[sin a hl{t -tc)+ cos a hl{t-tc)] , (3.1) 

where a depends on the orbital phase and orientation 
of the binary system, tc is the time (at the detector) 
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when the binary reaches its inner-most stable circular 
orbit, and hi ^(t — tc) are the two polarizations of the 
gravitational waveform produced by an inspiralling bi- 
nary that is optimally oriented at a distance of 1 Mpc. 
An optimally-oriented binary system is one that lies on 
the detector's z-axis with its orbital plane parallel to the 
x-y plane, defined by the arms of the detector (i.e., di- 
rectly above or below the detector and orbiting on the 
plane of the sky). The effective distance Deff depends 
on the true distance r to the binary, its location in the 
sky relative to the detector, and its orientation. This de- 
pendence is, in part, caused by the non-uniform detector 
response over the sky. If the source is not optimally ori- 
ented, then Doff > r. The binary inspiral waveform can 
thus be parameterized (for a single detector) in terms 
of the component masses, the effective distance, and the 
signal phase. 



The rate at which neutron star binaries coalesce in our 
Galaxy can be estimated using the observed sample of 
binary pulsars. (See, for example, Ref. |3|.) This rate 
estimate can be extrapolated to extra-galactic distances 
(following Phinney [sj) by assuming that the coalescence 
rate is proportional to the formation rate of massive stars 
and that the primordial binary population in our Galaxy 
is typical. Since the rate of massive star formation is pro- 
portional to blue-light (B-band) luminosity, the number 
of coalescences contributed by another galaxy is deter- 
mined by the ratio of its blue-light luminosity to that 
of the Milky Way. The sample population for our anal- 
ysis used spatial and mass distributions from a Milky 
Way population produced by the simulations of Ref. Q 
with the spatial distribution described in Ref. [s^. Ad- 
ditional sources from the Large and Small Magellanic 
Glouds, treated as points^ at their known distances and 
sky positions, were also added. The number of sources 
was proportional to the absolute blue-light luminosity of 
the LMC and SMC, with correction factors applied to 
account for reddening and the lower mctallicity of these 
objects. The latter leads to lower neutron star forma- 
tion rates primarily due to weaker stellar winds, which 
in turn favor the formation of more massive compact ob- 
jects. With these corrections, the event rates from the 
Large and Small Magellanic Clouds are taken to be 11% 
and 2% of the Milky Way rate. We note that this popu- 
lation model may not be exactly accurate, but is repre- 
sentative of the current understanding of binary neutron 
star formation. 
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FIG. 2: Summary of detector status and sensitivity to the 
population of neutron stars described in Sec. lIIII as a function 
of sidereal time. For a given sidereal time, the upper panel 
shows the number of days during the run when at least one of 
the interferometers (HI or LI) was collecting scientific data. 
For reference, the vertical dotted line indicates 05:00 UTC 
(corresponding to midnight at Livingston) on September 01, 
2002. The lower panel shows the effective distance as mea- 
sured in Livingston [and defined by Eq. (13.111 1 to 10%, 50%, 
and 90% of the binary neutron star population described in 
Sec. mil The horizontal dashed lines show the average dis- 
tance at which an inspiral of 2 x 1.4M0 neutron stars, in the 
optimal direction and orientation with respect to each detec- 
tor, would produce a signal-to-noise ratio of 8, i.e. 176 kpc 
for LI and 46 kpc for HI. 



IV. TEMPLATE BASED TRIGGER 
GENERATION 



The data stream from each detector was searched 
for inspiral waveforms using matched filtering, i.e., by 
evaluating the correlation (with a frequency-dependent 
weighting to suppress noise) between the data and a tem- 
plate waveform for all possible coalescence times. We use 
templates for non-spinning binaries, so each waveform is 
identified by a mass pair / = (777.1,777.2), a phase a and a 
distance D^s as described above. The gravitational wave 
signals also obey the approximate relationship 



^ The angular diameters of the Large and Small Magellanic Clouds 
are 7 and 4 degrees, respectively. These are comparable to the 
best angular resolution that can be achieved in our analysis us- 
ing time of arrival information from two LIGO detectors to de- 
termine sky position information. The resolved variations of in- 
strumental response across the Magellanic Clouds is negligible in 
our analysis. 



Kit) = -^Hif) , (4.1) 

where / > and the Fourier transform q{f) is defined by 



9(/) 



^2-»/*g(i) dt 



(4.2) 
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We exploit the symmetry (|4.1() , which is exact within the 
stationary-phase approximation used in this analysis,^ 
to reduce computational overhead in searching over the 
phase a. If the detector's calibrated strain data is s(t) = 
n(t) + where is the instrumental strain noise 
and h[£) is a gravitational wave signal (if present), then 
the matched filter output for given masses / ~ {mi, 1112) 
is the complex time series 

z{t) = x(i) + iyit) 4 ^'^PJ//^ e^^'^'df (4.3) 

where S„{f) is the one-sided strain noise power spectral 
density. In this expression, x{t) is the matched filter 
response to the a = waveform ft,^ while y{t) is the 
matched filter response to the a = Tr/2 waveform hi. 
Matched filtering theory [s^ provides a simple way to 
search over the phase a: construct the signal-to-noise 
ratio (SNR) of the matched filter output, 

Pit) ^ M, (4.4) 
a 

where 

is the variance of the matched filter output due to de- 
tector noise. For stationary and Gaussian noise, p is the 
optimal detection statistic for a single detector. 

The waveform (|3.1|1 depends on the masses of the two 
companions, so a bank of templates that covers the ex- 
pected range of neutron star masses must be used |36j . 
We adopted a template bank that covers the mass range 
1-3 Mq for each companion. The discrete bank was de- 
signed to cause less than 3% loss in SNR due to param- 
eter mismatches between any waveform and the nearest 
template in the bank. The layout of the template bank 
depends on the noise power spectral density of the instru- 
ment. A single template bank was used in this analysis: 
banks were first generated for each instrument and the 
bank with the most templates (in this case, the one gen- 
erated for LI) was used. We checked that the resulting 
2110 templates covered the mass range with < 2% loss of 
SNR for LI and < 7% loss for HI. Waveforms with total 
mass below A.OMq incurred < 3% loss of SNR in both 
instruments. Using a single template bank allows easier 
comparison of inspiral candidates in the coincidence step 
of our analysis. 

To reject transient noise artifacts that may excite a 
matched filter, but do not accumulate SNR as a chirp 
signal would, we employed an additional time-frequency 
veto in which the contribution to the filter output z{t) 



^ The stationary-phase approximation to the Fourier transform of 
inspiral template waveforms was shown t o be sufficiently accurate 
for gravitational- wave detection in Ref. |34| . 



from p frequency sub-bands is compared to the expected 
contribution for the templates [ill l37| . The frequency 
sub-bands were chosen so that the expected chirp would 
produce an equal contribution to both the real and imagi- 
nary components of the filter output from each sub-band. 
The chirp for each sub-band is filtered to produce the p 
complex-quantities ziit) and the statistic is constructed 
as 

X'it)^^j2\'^it)-zit)/p\'. (4.6) 

In the presence of Gaussian noise alone, is chi-squared 
distributed with v = 2p — 2 degrees of freedom. In this 
analysis, we did not optimize over different values of p, 
but chose p ~ 8 which worked well. 

If a putative signal h(t) has masses which do not ex- 
actly match any template in the bank, then has a 
non-central chi-squared distribution with 2p — 2 degrees 
of freedom and a non-central parameter A = 2p^e^, where 
p is the SNR for the signal and e is the fractional loss of 
SNR due to parameter mismatch. While it is possible 
to construct constant confidence thresholds on the non- 
central chi-squared distribution for various signals, in this 
analysis we simply require 

< 5(p + 0.03p2) (4.7) 

for any inspiral event, where p = 8 as described above. 
We refer to this cut as the x^-veto. Since the detector 
noise was not Gaussian, the threshold was selected based 
on performance in the playground data set described in 
Sec.lvland not using the exact result for the non-central 
chi-squared distribution. 

We identify possible inspirals in a single detector (HI 
or LI) by finding maxima of p{t) above a certain thresh- 
old (chosen to be p* — 6.5 in this analysis), subject to the 
X^-veto constraint of Eq. H4.7|l . and separated in time by 
at least the length of the template. Each such maximum 
is considered a trigger; the inferred coalescence time, p, 
and values are cataloged in a database along with 
the template parameters and effective distance (in Mpc) , 
-Dcff = a /p. 

Times when each interferometer was in stable opera- 
tion were identified as science mode epochs. These sci- 
ence mode epochs were analyzed in blocks of 256 sec- 
onds overlapped by 32 seconds as shown in Fig. 13 If 
there was not enough data at the end of a science mode 
epoch to take a 256 second block for analysis, the ex- 
tra data was dropped from the analysis. Each 256 sec- 
ond block was read by the LIGO Data Analysis System 
(ldas) [111 , which down-sampled it from 16 kHz to 4 kHz. 
The power spectrum of the data was estimated for each 
block by dividing it into four 64 second segments and 
taking the mean power spectrum of these four segments. 
The matched filter given in Eq. (|4.3|l was implemented on 
64 second data segments using routines in the LSC Algo- 
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Science Mode Epoch time 




FIG. 3: Times when an interferometer was in stable opera- 
tion were identified as science mode epochs indicated by the 
thick black lines at the top of the figure. These science mode 
epochs were analyzed in blocks of 256 seconds overlapped by 
32 seconds (indicated in white). If there was not enough data 
at the end of a science mode epoch to take a 256 second block 
for analysis, the extra data was dropped from the analysis. 
Each of these blocks were further divided into 7 overlapping 
segments of 64 seconds which were then searched for inspiral 
signals. The overlaps are needed to avoid contamination in 
the correlation used to compute the SNR. 

rithm Library (lal) [s^.^ In order to avoid end-effects 
in performing the correlation described by Eq. H4.3|l . we 
modified 1/S'„(/) so that its inverse Fourier transform 
had a maximum duration of ±16 seconds. The first 
and last 16 seconds of each filtered 64 second segment 
were ignored as corrupted by the end-effects of the filter. 
The 64 second segments were overlapped by 32 seconds — 
thus forming 7 overlapping segments in each 256 second 
block — so that no data was lost within each block. Since 
the blocks were also overlapped by 32 seconds, only the 
first 16 seconds of data from the first block and the last 
16 seconds of data from the last block were lost from 
each science-mode epoch. These effects combined result 
in the loss of 14 hours of data from each of the LI and 
HI interferometers. 

When the interferometers at Hanford and Livingston 
were in stable operation, wc checked for coincident signals 
to improve confidence in a detection. Since the Hanford 
and Livingston detectors are approximately co-aligned, 
they should observe essentially the same gravitational- 
wave signal.^ Ignoring mis-alignment and assuming the 
instrumental noise is Gaussian and uncorrelated, the op- 
timal detection statistic can be written as 

Pcohoront(^) = max (4.8) 

where and ZHi(i) ai'c the complex matched filter 

outputs from the LI and HI detectors, af^^ and are 



^ The analysis was performed on the MEDUSA comput- 
ing cluster at the University of Wisconsin-Milwaukee 
http : //www. Isc-group .phys .uwm. edu/beowulf /medusa 

* The two LIGO interferometers HI and LI are not exactly aligned 
due to the curvature of the earth. The effect of this curvature 
is to introduce small differences in response of each instrument 
to a real gravitational wave. We have ignored this effect at the 
present time, but plan to include it in future analyses. 



the variances of these matched filter outputs for the two 
detectors, r is the difference in the arrival time of the 
signal between the two detectors, and the maximization 
is performed over all possible values of r up to the light- 
travel time between the two detectors (±10 ms) pOll4l|. 
This statistic uses the same template in each instrument 
and assumes that the time of arrival is consistent with the 
light travel time between the instruments. Since a [Eq. 
(|4.5() ] depends on the inverse power spectral density, a 
large value indicates good sensitivity. If, for example, LI 
is considerably more sensitive than HI (as it was during 
SI), then (Tli ^ crm- Thus, one has |zli| ^ l-^Hil both 
during typical operation and when a signal is present, 
and a good approximation to the coherent statistic is 

Pcohcrcnt - kLlP/c^Ll = Phi ■ (4-9) 

Since LI was much more sensitive than HI during the 
SI run, ^coherent for an event seen while both detectors 
were operating is well approximated by the p value for 
LI alone; when only HI was operating, Pcohorent reduces 
to the p value for HI since the contributions from LI 
vanish. We also note that a binary inspiral signal would 
have pli ^ '^Pm, so a genuine signal would not produce 
a trigger in HI unless it appears in LI with very high 
SNR (greater than ~ 26). 



V. DATA QUALITY CRITERIA AND VETOES 

The performance of the LIGO interferometers varied 
significantly during the SI run on both long and short 
time scales. We omitted intervals of data from a given 
interferometer if it was not properly calibrated or if it 
had an unusually high level of noise, as described below. 
We also were able to veto some individual triggers which 
had a clear instrumental origin. To avoid statistical bias, 
the specific veto criteria were decided based on studies 
of a playground data set comprising roughly 10% of the 
data collected when all three interferometers were oper- 
ating. This data was excluded from calculation of the 
final analysis results. 



A. Instrumental calibration 

As mentioned in Section m the time variation of the 
interferometer response was tracked by continuously in- 
jecting sinusoidal signals with known amplitudes. The 
calibration was updated once per minute, and the anal- 
ysis of each 256-second block of data used the first avail- 
able calibration update within the block. There were 
periods of time when the sinusoidal injections were ab- 
sent, however, and the calibration could not be updated. 
Blocks of data in which such a calibration drop-out oc- 
curred were not analyzed. There were also some periods 
of time when HI calibration information was present but 
was deemed unreliable; these periods also were omitted 
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from the analysis. In total, 17 hours of HI data and 8 
hours of LI data were omitted from the analysis because 
of missing or unreliable calibration data. 

B. Noise level 

The noise in the gravitational-wave channel of each 
interferometer was sensitive to optical alignment, servo 
control settings, and environmental conditions. During 
most of the run, the noise level varied by less than a 
factor of two; however, there were a number of times 
when the noise level was significantly higher. We chose 
to omit these periods when the noise was particularly 
high. The specific criteria were developed by the work- 
ing group searching for gravitational-wave bursts and 
adopted for the inspiral analysis as well. Each inter- 
ferometer's performance was tracked by calculating the 
band-limited root-mean-square noise (BLRMS) in four 
frequency bands {Bi, B2, B3, ={320-400 Hz, 400- 
600 Hz, 600-1600 Hz, 1600-3000 Hz}. For each band, 
the noise power Pi{t) was calculated every 1/8 seconds, 
then averaged over 360-second time intervals and com- 
pared to the mean value Pi for all science-mode data 
collected. Based on empirical studies of correlations be- 
tween the power in each band and non-stationarity of 
the noise, we decided to eliminate any contiguous epoch 
of science data if there was any 360-second interval dur- 
ing the epoch for which Pi > 10 Pi or Pj > 3Pj for 
j = 2,3,4. This BLRMS cut removed 13 hours (8%) of 
the LI data and 43 hours (18%) of the HI data. 

Since the BLRMS cut uses the noise in the 
gravitational-wave channel to identify times when data 
quality is suspect, a sufficiently strong inspiral signal 
could potentially cause the veto to be invoked. Based 
on the known amplitude response of the instruments, 
we determined that a binary neutron star inspiral signal 
would be vetoed in this way only if it were closer than 
^ 300 pc, corresponding to a SNR of 4.7 x 10"^ in LI. 
By way of confirmation, we also computed Pi for periods 
when large-amplitude simulated inspiral waveforms were 
injected into the interferometers. The observed safety 
margin was consistent with the model calculations. Since 
^ 1% of the target population is within 300 pc of Earth, 
the systematic effects of the BLRMS cut on our search 
were negligible. 



C. Instrumental vetoes 

The data quality cuts described above addressed per- 
formance variations over long time scales. Each of the in- 
terferometers also exhibited non-stationary behavior on 
short time scales, with occasional glitches and/or brief 
periods of elevated broadband noise in the gravitational- 
wave channel. Because the matched filtering technique 
used in this analysis assumed the noise spectrum to be 
stationary over periods of several minutes, these tran- 



sients tended to excite the inspiral filter bank in such a 
way as to be recorded as triggers with fairly large SNR, 
even though they did not closely resemble the waveform 
of an inspiral. The veto [Eq. (|4.7|l ] eliminated many of 
these triggers, but some remained, appearing as a high- 
side tail in the SNR distribution of inspiral triggers found 
in the playground data set. 

We attempted to identify environmental or instrumen- 
tal origins for these high-SNR triggers by checking for 
coincident transients in the many auxiliary data chan- 
nels which were recorded along with the gravitational- 
wave channel. These included environmental monitoring 
sensors (seismometers, accelerometers, magnetometers, 
etc.) as well as various signals related to the operation 
of the interferometers. We evaluated several transient- 
detection algorithms, eventually choosing a simple one 
which applies a high-pass filter to the data and records 
excursions from zero which exceed a given size threshold. 
We developed an automated procedure to veto any inspi- 
ral trigger within a given time window around auxiliary- 
channel glitches found by this algorithm. For each of 
several promising auxiliary channels, the excursion size 
threshold and time window were tuned using the play- 
ground data set to maximize the number of triggers ve- 
toed without introducing undue dead-time. The results 
of these studies for each interferometer are summarized 
below. 

The HI detector experienced distinct glitches in the 
gravitational-wave channel at a rate of about 4 per hour. 
Although no external environmental cause was identified, 
nearly all of these glitches were clearly visible in an aux- 
iliary channel derived from a photo-diode at the interfer- 
ometer's reflected port. This channel is sensitive to the 
average arm length and is used to control the frequency 
of the laser light. We vetoed inspiral triggers within a 
±1 second window on either side of glitches found in this 
auxiliary channel; this veto condition introduced a dead- 
time of 0.2%. Based on the detector design, a real grav- 
itational wave would not be expected to appear with a 
significant amplitude in this auxiliary channel; we verified 
this experimentally by injecting simulated inspiral wave- 
forms into the interferometer arm length control servo 
(changing the arm lengths using electromagnetic actua- 
tion to push the suspended mirrors) and observing the 
signal strength in this and other auxiliary channels. 

High-SNR inspiral triggers in the LI detector were 
strongly correlated with transients in an auxiliary chan- 
nel derived from the photo-diode at the interferometer's 
antisymmetric port, nominally orthogonal in demodu- 
lation phase relative to the gravitational-wave channel. 
This auxiliary channel was not used to control any de- 
gree of freedom in the interferometer; it was sensitive 
to imbalance in the modulation sidebands and to align- 
ment fluctuations. This suggested its use as a veto chan- 
nel. Unfortunately, simulated inspiral waveforms injected 
into the arm length control servo appeared with non- 
negligible amplitude in this auxiliary channel. We sus- 
pect this was an artifact of injecting a large signal with 
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imperfectly balanced mirror actuators, introducing an os- 
cillatory misalignment. To be safe, however, we chose not 
to veto based on this channel. No other auxiliary chan- 
nel offered an efficient veto, so no instrumental veto was 
applied for LI. 



VI. ANALYSIS PIPELINE AND TUNING 

The detection of a gravitational-wave inspiral signal 
in the SI data would (at the least) require triggers in 
both LI and HI with consistent arrival times (separated 
by less than the light travel time between the detectors) 
and waveform parameters. Such a temporal coincidence 
requirement has the advantage of greatly reducing the 
background rate due to spurious triggers in the individual 
detectors. It limits the volume of space searched to that 
which can be seen by the less sensitive detector, however, 
and it limits the observation time to the periods of simul- 
taneous operation. Because the LI detector was much 
more sensitive than HI during the SI run, and because 
they operated simultaneously less than 30% of the time, 
we developed a more sophisticated (upper-limit) analysis 
pipeline which makes use of triggers from the individual 
detectors when a coincidence test is not possible. Studies 
of the playground data set indicated that the additional 
background rate introduced by this choice should not off- 
set the improvement in event rate limit that comes from 
increased observation time. Of course, event candidates 
identified during non-coincident observation times could 
not lead to an unambiguous detection of gravitational 
waves. 

Our analysis pipeline is summarized in Fig. ^ We 
follow five steps to produce a list of non-vetoed event 
candidates which represent the background due to detec- 
tor noise (plus any gravitational-wave signals, if present) 
during periods of nominal operation. (1) Analyze the 
gravitational-wave channel data from each detector us- 
ing matched filtering as described above. When p > 6.5 
in an individual detector, apply the veto to eliminate 
spurious excitations of the templates. Store information 
about the surviving triggers in a database. (2) Apply the 
BLRMS cut to reject triggers in periods with unusually 
high noise, and apply a veto to eliminate HI triggers with 
a clear instrumental origin. (3) When both interferome- 
ters are operating, require coincident triggers only if the 
effective distance measured by the LI detector is closer 
than a cut-off distance D*. (The selection of D* and the 
coincidence criteria is described below.) In this case, the 
SNR for the event candidate is taken to be the LI SNR 
in accordance with the discussion around Eq. (|4.9() . If 
an LI trigger has Dcs > D*, keep the trigger regardless 
of whether it was also detected by HI. (4) During times 
when only one interferometer is operating, keep any trig- 
ger that passes the cuts in the second step. (5) Finally, 
maximize all surviving triggers over time and over the 
template bank. The timing resolution of inspiral signals 
is ^ 1 ms once coincidence of template mass parame- 
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FIG. 4: The inspiral analysis pipeline used to determine the 
reported upper limit. "HI Only", "HI & LI", and "LI Only" 
indicate which interferometer(s) was/were operating when a 
trigger was recorded. This method of recording candidate 
events even when coincidence is not available allows a tighter 
bound to be placed on the rate of binary neutron star inspirals 
by providing more observation time and allowing for the much 
greater sensitivity of LI than HI. 



ters in both instruments is enforced. When coincidence 
is unavailable, background noise can trigger many tem- 
plates at significantly different times. Since the impulse 
response of the matched filter is ~ 16 seconds [because 
the template is effectively convolved with the frequency 
dependent weighting l/5'„(/) when computing the SNR 
in Eg. I4.3j . we maximize over all triggers in a 16 second 
window and over the entire template bank to produce 
the final list of candidate events. The post-processing 
analysis described by steps (2)-(5) was performed using 
software in the package lalapps |39| . 

We characterized our analysis pipeline using a Monte 
Carlo method in which we re-analyzed the data with sim- 
ulated inspiral signals injected into the time series. The 
re-analysis used exactly the same pipeline as the origi- 
nal analysis and the simulated signals were drawn from 
the population described in Sec. 11111 The efficiency of 
the pipeline is the fraction of this population that could 
be detected. To avoid statistical bias, we used only the 
playground data set described in Sec. when deciding 
aspects of the pipehne. 

The coincident event selection criteria in step (3) were 
tuned by studying the fractional loss of efficiency of the 
pipeline. A trigger from HI was considered coincident 
with a trigger from LI if the recorded coalescence times 
were within a time window At* = 0.011 s. This accounts 
for the light travel time between the two sites (which is 
0.010 s) plus statistical and systematic errors in the in- 
dividual measurements of coalescence time. The gross 
frequency evolution of an inspiral chirp signal is con- 
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TABLE I: Number of event candidates with Pcohcront > 8.0 
found via each of the pipeline paths shown in Fig. 01 The 
first two lines represent event candidates found while both 
interferometers were operating. No coincident events were 
detected in both interferometers; however, there were many 
event candidates found in LI with effective distances Deff > 
51 kpc, which would not be detectable in HI and thus are 
kept as event candidates. The last two lines represent event 
candidates found while only one interferometer was operating. 



Operating Detected in Number Max SNR 

LI and HI LI (Detr < 51 kpc) and HI 

LI and HI LI {D^b > 51 kpc) 418 15.6 

LI only LI 786 15.9 

HI only HI 274 12.0 



trolled by the chirp mass Ai = m^^^m^^^ {rai +777.2)"^^^. 
The difference of chirp mass = AIli — A^Hi for a 
pair of coincident (in time) triggers was required to sat- 
isfy \/S.M\/M'Li < 10"^ leading to ~ 1% fractional loss 
of efhciency for the playground data. Finally, we chose 
D* = 51 kpc, producing 10% fractional loss of effi- 
ciency for the playground data, in order to have a rea- 
sonable chance of detection in coincidence between the 
two sites. 



are summarized in Table 

No event candidates were found in coincidence by both 
detectors. If there had been one or more coincident 
event candidates, the background rate of accidental co- 
incidences could have been determined from the data by 
counting coincidences after shifting the HI trigger times 
relative to the LI trigger times by an amount greater 
than the light travel time between the sites. In fact, 
in the SI data, there were no triggers whatsoever in LI 
which were close enough (Z^off < 51 kpc) to have been 
seen in HI with pni > 6.5. 

For comparison. Table shows the number of events 
identified with Pcohorcnt > 8.0 by the same analysis 
pipeline upon processing the output of the Monte-Carlo 
simulation described in Sec. IVII A total of 5071 sim- 
ulated signals were overlaid on the SI data, of which 
619 were found in coincidence, demonstrating that the 
pipeline could correctly identify coincident event candi- 
dates within 51 kpc. Note that the counts of event can- 
didates in the other three paths of Table IhI include those 
in the underlying data, not associated with an injected 
signal. 

TABLE II: Results from the Monte-Carlo simulation given 
for comparison with the equivalent results of the search. 
Note that 619 simulated events were detected in coincidence, 
demonstrating that the pipeline was indeed capable of iden- 
tifying coincident event candidates. 



VII. RESULTS FROM SI DATA 

The non-playground data was analyzed using the 
pipeline described above. After the division of the data 
into 256-second blocks, the rejection of blocks without re- 
liable calibration, the additional loss of 16 seconds from 
the beginning of the first block and the end of the last 
block of a science-mode epoch, and the times during 
which a veto was active were discarded, a total of 236 
hours of non-playground data remained: 58 hours when 
both LI and HI were operating, 76 hours when only LI 
was operating, and 102 hours when only HI was operat- 
ing. 



A. Triggers and Event candidates 

The triggers from each interferometer satisfy 
/^coherent > 6.5 and the veto defined in Eq. H4.7|l . 
There were ~ 2 x 10® triggers from each detector 
before applying vetoes, checking for coincidence, and 
maximizing over templates and time with a 16 second 
window. The numbers of event candidates from each 
part of our pipeline with Pcohercnt > 8.0 in the SI data 



Operating Detected in Number Max SNR 



LI and HI 


LI (Oeff 


< 51 kpc) and HI 


619 


634.4 


LI and HI 


LI (Deff 


> 51 kpc) 


773 


46.5 


LI only 


LI 




2052 


460.2 


HI only 


HI 




1623 


221.9 



The ten event candidates with the largest SNR in the 
pipeline were all detected by LI and had SNR between 12 
and 16 and per degree of freedom between 2.2 and 4.9. 
Details of the five largest events are given in Table ITTTl 
Four of these events have values close to the thresh- 
old in Eq. H4.7|) : the exception is the candidate which 
occurred at 13:06:56.731 UTC on 2002/09/02. FigureEl 
(left panels) shows the signal-to-noise and time series 
for the candidate with the largest SNR, which occurred 
at 00:38:33.557 UTC on 2002/09/02. A simulated inspi- 
ral signal with comparable SNR is shown in Fig. |S1 (right 
panels) to demonstrate the qualitative differences in the 
time series. Unlike the simulated signal, the SNR of the 
event candidate is consistently high across the duration 
of the event, with the value of the veto varying sig- 



^ Since our pipeline with Pcoherent > 6.5 identifies a high num- 
ber of candidate events (close to the maximum number possible 
for our pipeline choices), we show only candidate events with 
Pcoherent > 8-0 in Tables HI and ITTl 
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TABLE III: The five candidates witli tiie largest SNR wiiicii remain at tlie end of tlie pipeline. This table indicates the time 
they registered in the detectors, the SNR, the value of per degree of freedom, the effective distance to an astrophysical event 
with the same parameters, and the binary component masses of the best matching template. 



Date UTC GPS Time Operating Detected in SNR x^/d-O-f- -Pcff (kpc) mi (Mq) m2 (Mq) 



2002/09/02 


00:38:33.557 


714962326.557 


LI only 


LI 


15.9 


4.3 


95.0 


1.31 


1.07 


2002/09/08 


12:31:38.282 


715523511.282 


LI and HI 


LI (L»efi > 51 kpc) 


15.6 


4.1 


68.4 


1.95 


0.92 


2002/08/25 


13:33:31.000 


714317624.000 


LI only 


LI 


15.3 


4.9 


100.7 


3.28 


1.16 


2002/08/25 


13:29:24.250 


714317377.250 


LI only 


LI 


14.9 


4.6 


88.7 


1.99 


1.99 


2002/09/02 


13:06:56.731 


715007229.731 


LI only 


LI 


13.7 


2.2 


96.3 


1.38 


1.38 



nificantly and dropping below the threshold right at the 
time of maximum SNR. 

Further scrutiny of the five largest SNR events revealed 
some instrumental problems. The event at 00:38:33.557 
UTC on 2002 /09 /02 coincides in time with saturation of 
the photo-diode at the antisymmetric port. This satu- 
ration, which started a second before the recorded coa- 
lescence time for the candidate event and lasted several 
seconds, was likely due to an instrumental misalignment. 
The misalignment is indicated by a five-fold increase in 
the power at the dark port of the interferometer, start- 
ing three seconds before the coalescence time and lasting 
six seconds. This event would have been vetoed by the 
auxiliary-channel veto condition we considered for LI but 
decided not to use (as discussed in Sec. IV C|) . The event 
recorded at 13:06:56.731 UTC on 2002/09/02 occurred 
when the interferometer was kept functioning during the 
most severe seismic conditions for SI data. Another event 
candidate, with SNR 13.0, occurred just 98 seconds later. 
The interferometer was rarely locked with seismic noise 
this high, and was probably experiencing up-conversion 
of low- frequency seismic noise into the gravitational- wave 
band through coupling with mechanical resonances and 
power line harmonics. 

Event candidates detected in just one interferometer 
cannot be taken to be real gravitational wave inspirals 
with any confidence, since we do not understand the dis- 
tribution of background. However, we can still place an 
upper limit on the rate of inspirals. Despite being able 
to find a posteriori reasons to justify eliminating some of 
the largest SNR event candidates as instrumental effects, 
we chose to keep them as event candidates for purposes 
of calculating the upper limit. 



B. Upper Limit Analysis 

To determine an upper limit on the rate of binary neu- 
tron star inspirals, we compare the observed distribution 
of events as a function of /Ocohcrcnt to the expected back- 
ground plus the population of interest. The comparison 
is made based on criteria established in advance of the 
analysis. Typically, one might choose an SNR thresh- 
old p* based on the rate and distribution of background 



events and compare the number of observed events with 
p > p* to the expected background. Unfortunately, we 
have no model for the background events in each of the 
interferometers; this is problematic because we chose to 
include event candidates found in only one interferome- 
ter to increase the visible distance and observation time. 
Rather than choosing a fixed value for p* , we adopt an 
approach in which p* is determined by the data. Specif- 
ically, we set p* equal to the largest SNR observed in 
the data and calculate the efficiency of the pipeline ac- 
cordingly. Since no events are observed with p > p*, we 
calculate an upper limit on the event rate for the mod- 
eled population assuming the probability of a background 
event above this SNR is negligibly small. This approach 
has the advantage of dealing with the lack of a model for 
the background events in a controlled manner. 

If the population of sources produces Poisson- 
distributed events with a rate TZ, the efficiency e(p*) is 
also the probability that any given binary neutron star 
inspiral in the target population would have SNR greater 
than p*. Then the probability of observing an inspiral 
signal with p > p*, given some rate TZ and some obser- 
vation time T, is 

P{p> p*-n) = l-e-'^'^<p'\ (7.1) 

A frequentist upper limit with 90% confidence on the 
value of TZ is determined by solving P{p > Pmax; '^90%) = 
0.9 for 7?.9o% where pmax is the largest SNR event ob- 
served in the SI data. The result can be written in closed 
form as 

2.303 

7^90% = (7.2) 

where emax = e(Pmax) and T is the observation time. 
For TZ > TZgQ%, there is more than 90% probability that 
at least one true inspiral event would be observed with 
SNR greater than Pmax- This limit is conservative since 
the non-zero probability that a background event could 
have SNR greater than Pmax has been neglected. 

It is useful to express the limit as a rate per Milky- Way 
Equivalent Galaxy (MWEG) for easy comparison with 
theoretical predictions and other observational results. 
The effective number of Milky Way equivalent galaxies 
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FIG. 5: Left Panels: The largest SNR candidate event seen during our search of the LIGO data. This candidate event 
occurred at a time when only LI was in stable operation. The top panel shows the signal-to-noise time series, p{t). Notice 
that p{t) > 6.5 many times in a ~ 5 second interval around the candidate event. The center panel shows x^/iP + 0.03p^) as 
a function of time; notice /(P + 0.03p^) > 5 for ~ 5 seconds around the candidate event, but drops below this threshold 
right at the time of maximum p. The inset shows this more clearly for ±0.1 second around the event where the threshold is 
indicated by a dot-dashed horizontal line. The bottom panel shows the time series for this candidate event after applying a 
high-pass filter with a knee frequency of 200 Hz. Notice the bursting behavior which does not look like an inspiral chirp signal. 
Right Panels: A simulated injection into the LI data. This example was chosen for comparison with the largest SNR event 
shown in the left panels since it similar in mass parameters, detected signal to noise and ■ The instrument was behaving 
well at the time around the simulated injection. The top panel shows that p{t) < 6.5 except in close proximity to the signal 
detection time. The center panel shows X^/(P + 0.03p^) as a function of time. Notice that it is much closer to threshold at 
all times around the simulated injection; this contrasts dramatically with the case of the candidate event shown in the left 
panels. The inset shows this more clearly for ±0.1 seconds around the injection. The bottom panel shows the time series for 
this simulated injection after applying a high-pass filter with a knee frequency of 200 Hz. The inspiral chirp signal is not visible 
in the noisy detector output. 



to which the search was sensitive is 

Ng = (^:^) (7.3) 

where Lq = 9 x IQ^Lq is the effective blue-light luminos- 
ity of the Milky Way and Lpop is the effective blue-light 
luminosity of the population. The rate limit can be writ- 
ten as 

n,o% = 2.303 X (^^^ (^^^ y-i MWEG-i . (7.4) 

During the T — 236 h = 0.027 y of data used in our anal- 
ysis, the largest observed SNR was Pmax = 15.9. The 
detection efficiency was computed using a Monte Carlo 
simulation in which we re-analyzed the data with sim- 
ulated inspiral signals, drawn from the population de- 
scribed in Sec. IIIII injected into the time scries. The 



efficiency e(p*), shown in Fig. |H{b), is the fraction of the 
5071 simulated signals which were detected with p > p*. 
The efficiency at p* = 15.9 is Smax = 0.53. Folding this 
together with Lpop = 1.13-Lg; the nominal value of A^g 
is 0.60; however, this is subject to some uncertainties, to 
be discussed in the next section. As a function of the 
true value of Nq, the rate limit is 

7^9o% = 1.4 X 102 (^^^ y^i MWEG-i . (7.5) 

It is interesting to compare our result with a direct 
estimate based on average sensitivity of the instruments 
(as shown in Fig. [21), properties of the population, and the 
observation times used in this analysis. At SNR 15.9, LI 
was sensitive to 80% of the sources and HI was sensitive 
to 35% of sources in our model population. Out of 236 
hours, LI was the best detector for 134 hours and HI for 
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FIG. 6: Panel (a) shows the number of events in the data with 
SNR > p* as a function of p* . The largest event has SNR = 
15.9. Panel (b) shows the detection efficiency e(/9*) for sources 
in the target population (Milky Way and Magellanic Clouds) 
as a function of p* . The dashed lines indicate boundaries of 
our estimated systematic errors on the efficiency. 

102 hours. The expected efficiency is then 

e(15.9) = (102 X 0.35 + 134 x 0.80)/236 = 0.6 . (7.6) 

The measured efficiency is e(15.9) — 0.53, but the veto 
and coincidence requirements both introduce some loss; 
the expectation based on playground data was w 0.06 x 
58/236 = 0.015 decrease in efficiency from coincidence 
and a loss of about « 0.06 from the . The actual loss 
from coincidence is « 0.02 as measured on the full data 
set. Consequently, the measured efficiency and hence the 
upper limit agree well with expectations. 

VIII. ERROR ANALYSIS 

The interpretation of this search for gravitational 
waves from binary neutron star inspiral suffers from a 
number of systematic effects which could modify the up- 
per limit. We classify these effects into three different 
types: (i) uncertainties in the population model and the- 
oretical expectations about the sources; (ii) uncertainties 
in the instrumental calibration; (iii) deficiencies of the 
analysis pipeline. Each one can have a direct effect on 
the efficiency of the search to detect gravitational waves 
from the target population as it exists in nature. 

A. Uncertainties in population model 

Uncertainties in the population model used for the 
Monte-Carlo simulations may lead to differences between 
the inferred rate and the rate in the universe. Since the 
effective blue-light luminosity ipop is normalized to our 



Galaxy, variations arise from the relative contributions 
of other galaxies in the population. These contributions 
depend on the estimated distances to the galaxies, esti- 
mated reddening, and corrections for metallicity (lower 
values tend to produce higher mass binaries), among 
other things. Since the Magellanic Clouds contribute 
only 13% of the blue light luminosity in this anal- 
ysis, a conservative estimate of the uncertainties gives 
Lpop = 1.13 ±0.06. 

The spatial distribution of the sources can also intro- 
duce significant uncertainties. Typically, the distances to 
nearby galaxies are only known to about 10% accuracy. 
Uncertainties in distances to galaxies near the limit of 
detector sensitivity are most relevant. As the detector 
sensitivity improves, more galaxies will be in this cat- 
egory, so it may become a major source of systematic 
uncertainty. It is not important for the current analy- 
sis, since the detectors were sensitive to the majority of 
sources in the Milky Way and Magellanic Clouds. 

The effects of spin were ignored both in the popula- 
tion and in the waveforms used to detect inspiral sig- 
nals. Apostolatos has performed the most complete 
analysis of the effects of spin on detection of waves from 
neutron star inspiral. His investigations suggest that less 
than 10% of all possible spin orientations cause more than 
~ 5% reduction in SNR for binary neutron star systems. 
There is insufficient information about the distribution 
of binary spin orientations to quantitatively estimate the 
systematic effect, but it seems certain that the fraction 
of the population with spin configurations which would 
interfere with their detection is negligible. 

Different models for NS-NS formation can lead to small 
variations in the tails of the NS mass distribution |^, 
but the bulk of the distribution always remains strongly 
peaked around observed NS masses |43|. Since the de- 
tection efficiency depends most sensitively on the bulk 
properties of the mass distribution, the expected varia- 
tions are negligible compared to other systematic effects 
discussed in this section. 



B. Uncertainties in the instrumental response 

The instrument response i?(/) was constructed for ev- 
ery minute of data during SI from a reference sens- 
ing function C(/), a reference open loop gain function 
G(/), and a parameter ait) representing varying opti- 
cal gain |23|. The parameter a was reconstructed using 
the observed amplitudes of the calibration lines described 
in Section^ If an inspiral signal is present in the data, 
systematic errors in the calibration can cause a mismatch 
between the template and the signal. For simulated in- 
jections, the SNR differs from the SNR that would be 
recorded for a signal from a real inspiral event at the 
same distance as the injection. The effect is linear in 
amplitude errors causing either an upward or downward 
shift in SNR, but quadratic in phase errors causing an 
over-estimation of sensitivity. This effect is captured by 
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shifting the cfhcicncy curve in Fig. horizontally by the 
appropriate amount. 

A careful evaluation of uncertainties in the SI calibra- 
tion j2fll | has shown that amplitude errors are primarily 
due to statistical fluctuations in the measurement pro- 
cedure, while phase errors are mostly systematic and 
are greater at higher frequencies. Combining statistical 
and systematic errors in quadrature, the amplitude errors 
lead to ~ 18% errors in SNR in LI and ~ 8% errors in 
HI. The phase errors lead to overestimation of the SNR 
by 2% in LI and 4% in HI. Combining amplitude 
and phase errors in quadrature and taking the larger LI 
values as representative, we find ^ 18% errors in SNR 
of Monte-Carlo injections which translates to fractional 
errors in efficiency ~ +14%/-10%, i.e. e,„ax = 0.53t°a"J^. 

To verify the data analysis methods, a few special stud- 
ies were done in which simulated inspiral waveforms were 
injected into the interferometer hardware using the mir- 
ror actuators. We then used the analysis pipeline de- 
scribed above to recover the known mass and distance 
parameters of the injected signal. A side benefit of these 
injections is to build confidence in our understanding of 
calibration uncertainties. In order to simplify the anal- 
ysis pipeline, the template bank was reduced to a single 
template, a 1.4, 1.4M0 or a 4.O,1.4M0 inspiral, corre- 
sponding to the mass parameters of the injected signal. 
Unfortunately, the calibration signal was turned off dur- 
ing the injections, so we defined a set of possible response 
functions for this range, and studied the variation in the 
detected inspiral signal. This was possible because the 
parameter a has only a limited physical range. We found 
that the variation in the reconstructed signal to noise 
and effective distance was in agreement with our expec- 
tations. Since the parameter a has a known dependence 
on the interferometer alignment we were able to use aux- 
iliary channel information to estimate its value during 
the injections. For this value the detected coalescence 
time of the chirp was the same as the injected time to 
within 1/4096 seconds, i.e. one sample of filtered data, 
and the reconstructed distance and the injected distance 
agreed to within 12%, which is consistent with the errors 
quoted above. 



C. Uncertainties in the analysis pipeline 

Since we use matched filtering to search for gravita- 
tional waves from inspiralling binaries, differences be- 
tween the theoretical and the real waveforms could also 
adversely effect the results. These effects have been 
studied i n g reat detail for binary neutron star sys- 
tems HH El The results indicate 10% loss of 
SNR due to inaccurate modelling of the waveforms for 
binaries in the mass range of interest. This feeds into 
our result through our measurement of the efficiency. 
We may be over-estimating our sensitivity to real bi- 
nary inspiral signals; this would shift all points on the 
efficiency curve in Fig. Elto the left by 10%. This cor- 



responds to fractional errors ^ +0%/— 5% in efficiency. 

I.e. Cmax — —0.03' 

The effects of discreteness of the template placement, 
errors in the estimates of the power spectral density 
Sn{f) used in the matched filter in Eq. H4.3|) . and trends 
in the instrumental noise are all accounted for by the 
Monte-Carlo simulation. 



D. Combined uncertainties on A'^g and the rate 

The efficiency incurs fractional errors ~ -1-14%/— 10% 
from calibration uncertainties (Sec. IVHI Bp and ^ 
-1-0%/— 5% from inaccurate knowledge of the inspiral 
waveforms fSec. lVIII Cll . Combining these in quadrature 
yields total errors ~ -1-14%/— 11% in the efficiency emax- 
Adding these (not in quadrature) to the ±5% error for 
Lpop rSec lVlll At yields 

A^G = 0.60l°:ig . (8.1) 

To be conservative, we assume the downward excursion 
A^G = 0.60-0.10 0.50 when using Eq. lf73|l to derive an 
observational upper limit on the rate of binary neutron 
star coalescence: 

7^ < 1.7 X 102 y-i MWEG"^ . (8.2) 



IX. CONCLUSIONS AND FUTURE 
DIRECTIONS 

The first search for gravitational-wave signals from 
coalescing neutron stars in LIGO science data yielded 
no coincident event candidates. An observational upper 
limit 1.7 X 10^ y~^ MWEG~^ on the rate of neutron star 
inspirals was derived. This limit is better than previous 
direct limits by a factor of 26 . 

Over the next few years, the sensitivity of the LIGO in- 
terferometers will be dramatically improved, to the point 
where inspirals of double neutron stars (2 x 1.4AfQ) are 
expected to be detectable within an equivalent volume 
w (47r/3) X (21 Mpc)^ E3|. Due to the non-uniform re- 
sponse of the detectors, this implies that a neutron star 
inspiral could be detected out to a maximum distance 
« 46 Mpc if the binary is located directly above or be- 
low the detectors and orbits in the plane of the sky. The 
rate of coalescence of extra-galactic neutron star binaries 
is thought to be proportional to the rate of massive star 
formation which is, in turn, proportional to the blue light 
luminosity. (See, for example, Ref. 0.) Using current 
galaxy catalogs, it is estimated that Nq « 465 MWEG 
will be detectable by LIGO (using the three detectors 
combined to produce a network SNR > 8) If the 

coalescence rate of binary systems (in which each com- 
ponent has a mass in the range 1-3 Mq) were as high 
as ~ 5 X 10-4 y-i MWEG"^ 0,113, then the event rate 
detectable by LIGO would be A''g times higher providing 
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up to 1/4 events per year. In lieu of a detection, an up- 
per limit within the range of astrophysical expectations 
will constrain the binary neutron star population mod- 
els, and especially the population of electromagnetically- 
undetectablepulsars at the faint end of their luminosity 
function Is^. 

The methods used, and experience gained, on the 17- 
day SI data set will be enhanced and used in future 
searches for gravitational waves from coalescing compact 
binaries with LIGO data. We can expect improvements 
in the upper limits obtained with detectors of better sen- 
sitivity, but we can also draw lessons on the methods 
used from this first experience. For example, we expect 
to reduce the maximum SNR of non-gravitational wave 
signals by making better use of the knowledge of the in- 
strument status to find more efficient veto criteria. In 
our next search, wc will require coincidence from candi- 
dates from the two observatories to establish an event. 
This will allow us to to measure a background rate of ac- 
cidental coincident events, using techniques to find lower 
SNR triggers as needed in the least sensitive instrument 
(if there continue to be significant differences in sensitiv- 
ities). Eventually, wc would like to use coherent methods 
with all the different detectors in operation. Even though 
the errors in the upper limits obtained in this article do 
not compromise their significance, the same errors would 
affect more seriously the parameter identification of a de- 
tection, so we hope to improve on all aspects contributing 
to statistical and systematic errors. 

Future searches will also target neutron-star/black- 
holc and black- hole/black- hole binaries which produce 
more energy in gravitational waves and will be visible 
within a much greater volume of the Universe. It is pos- 
sible that several black-hole binaries could be detected 
by the initial LIGO interferometers |^, ^3 1 but there is 
considerable uncertainty in this event rate. An obser- 
vational upper limit would constrain population models 
and yield information about the formation mechanisms 
of black-hole binaries. The challenge of setting an upper 
limit on higher-mass binary systems is formidable: mas- 
sive binary systems (black-hole/black- hole) will exhibit 
highly relativistic effects (beyond the realm of the stan- 
dard post-Newtonian approximation) within the sensitiv- 
ity band of the instruments |29l . , whereas spin-orbit 
and spin-spin coupling in precessing binaries will be ex- 
tremely important in intermediate-mass systems of low 
mass ratio (neutron-star/black-hole) |3lL liM ISOl ISlL IS ^ . 



These effects will greatly expand the parameter space 
that needs to be searched, and will require the construc- 
tion of both accurate and computationally efficient 
waveforms. Efforts are already under way to construct 
detection template families |48l l53l Is^ in our search 
codes. The goal with these detection template fami- 
lies is to efficiently mimic all the known analytical mod- 
els of black-hole binary dynamics (such as the standard 
post-Newtonian models 1281 and their improved versions, 
namely ^ P-app roximants |55j and effective one-body tech- 
niques [2^ |5a 113,111, 113) and/or the effects of precession 
on waveforms emitted by binaries with spinning compact 
objects. Despite the challenges, a search for gravitational 
waves from black hole binaries is the highest priority for 
current research. 

Another class of systems is the sub-solar mass 
(0.2-1 A/0) binary black holes that might form a sizable 
portion of macroscopic halo objects (MACHOs) 113 • If 
such objects exist, then many of the challenges in detect- 
ing binaries with stellar-mass are alleviated: the orbits 
of these binaries will not be highly relativistic while the 
gravitational waves are emitted in the LIGO sensitivity 
band, and the spin effects can be handled easily. On 
the other hand, the smaller amplitude of the gravita- 
tional waves emitted by these sources limits the distance 
to which they can be seen. 
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